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We study the K* meson reduction in heavy ion collisions by focusing on the hadronic effects on 
the K* meson abundance. We evaluate the absorption cross sections of the K* and K meson by light 
mesons in the hadronic matter, and further investigate the variation in the meson abundances for 
both particles during the hadronic stage of heavy ion collisions. We show how the interplay between 
the interaction of the K* meson and kaon with light mesons in the hadronic medium determines the 
final yield difference of the statistical hadronization model to the experimental measurements. For 
the central Au+Au collision at ^/sNN = 200 GeV, we find that the K* jK yield ratio at chemical 
freeze-out decreases by 36% during the expansion of the hadronic matter, resulting in the final ratio 
comparable to STAR measurements of 0.23 ±0.05. 
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I. INTRODUCTION 

Relativistic heavy ion collision experiments have en¬ 
abled the production of a system of quantum chromody¬ 
namic (QCD) matter at extreme conditions under con¬ 
trolled conditions [IHI. Due to the huge energies avail¬ 
able in heavy ion collisions, it is expected that a possible 
phase transition predicted by Lattice QCD @ between a 
hadronic matter and a system of deconfined quarks and 
gluons takes place, and the quark-gluon plasma (QGP) 
at very high temperature is produced at the initial stage 
of the collision. As a result, large numbers of hadronic 
particles are produced during the quark-hadron phase 
transition at later stages of heavy ion collisions. 

These hadronic particles are believed to emerge at the 
transition point with the information of the matter. The 
statistical hadronization model has been quite success 
in explaining the measured production yields of hadrons 
with two parameters characterizing the chemical freeze- 
out point in heavy ion collisions; the phase transition 
temperature and the baryon chemical potential [M3- 

All particles produced at the freeze-out, however, are 
subject to further interactions with other hadrons in the 
hadronic matter, leading to possible deviations in the 
final yield of some hadrons from the statistical model 
prediction. In addition to the effects from hadronic in¬ 
teractions, the lifetime of hadrons as well as the lifetime 
of the hadronic matter itself plays an important role in 
changing the abundance of hadrons from the yield at the 
chemical freeze-out. 

The abundance of hadrons stable against strong de¬ 
cays is expected to be changed mostly by hadronic inter¬ 
actions while that of resonances will be affected by both 
their interactions with other hadrons and their strong de¬ 
cays when the lifetime of resonances is comparable to, or 
smaller than the lifespan of the hadronic stage in heavy 
ion collisions. Daughter particles of resonances are sub¬ 
ject to re-scatter as well in the hadronic medium, making 
the reconstruction of the resonances from an invariant 
mass analysis difficult. 


Studying the effects from the hadronic interactions on 
the abundance of resonances has been suggested as one 
way of confirming the scenario about a time delay be¬ 
tween the chemical and thermal freeze-out pm , since 
a sudden hadronization in heavy ion collisions would 
leave no time for resonances to decay in the hadronic 
medium. In particular, the K* meson has attracted lots 
of attentions as its short lifetime 4 fm/c is less than the 
presumed lifespan of the hadronic stage. 

The effects of hadronic interactions on the yield of the 
K* meson have been measured in heavy ion collisions us¬ 
ing K* jK yield ratios. Since the K meson is the ground 
state of the K* meson, having the same valence quarks 
with a different mass and relative orientation of its quark 
spins, the K*/K yield ratio is considered to be indepen¬ 
dent of the freeze-out conditions when hadronic interac¬ 
tions are neglected. It has been shown that K* /K yield 
ratios decrease with the increasing size of the system at 
the same energy PP. Compared to p±p collisions, 
K*/K yield ratios in Cu±Cu and Au±Au collisions are 
smaller, naively implying that K* and K mesons partic¬ 
ipate in re-scattering processes during the expansion of 
hadronic matter, and that the hadronic effects become 
larger as the size of the hadronic matter increases. 

The average transverse momentum of the K* meson 
measured in heavy ion collisions PP, which is higher 
than that of the K* meson in p±p collisions, also sup¬ 
ports the re-scattering scenario about K* mesons. K* 
mesons with low transverse momenta escape the hadronic 
stage later than K* mesons with higher transverse mo¬ 
menta, and thus suffer more re-scattering in the hadronic 
medium. As a result the measurement of the K*/K yield 
ratio 0.23 ± 0.05 P in Au±Au collision at ^/snn = 200 
GeV is smaller than the statistical model expectation 
0.33 ±0.01 at that collision P. However, as we will see, 
the measurement of the K* jK yield ratio inconsistent 
with the statistical model prediction not only confirms 
the hadronic effects on the yield of the K* meson but 
also provides information on the change in the proper¬ 
ties of the hadronic matter at freeze-out. 

With these measurements of the K*/K yield ratio in 
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mind, we study here the hadronic effects on the K* me¬ 
son by evaluating its absorption cross sections with pions, 
p, K, and K* mesons, and furthermore investigate vari¬ 
ations in the K* meson abundance during the hadronic 
stage of heavy ion collisions by solving a time evolution 
equation for the K* meson. After the K* meson is pro¬ 
duced at the chemical freeze-out, it interacts mostly with 
light hadrons during the expansion of the hadronic mat¬ 
ter. As a result, K* mesons can be absorbed by the 
co-moving light mesons, or additionally produced from 
scattering between them. Thus, evaluating the K* meson 
cross sections by light hadrons is necessary in estimating 
the hadronic effects on the K* meson abundance in heavy 
ion collisions. By comparing our results with the experi¬ 
mental observation in heavy ion collisions, we understand 
the discrepancy of the K* meson yield between the sta¬ 
tistical model and the experimental measurements. 

As has been stated, the scattering of the K* meson 
daughter particles such as kaons in the hadronic medium 
also contains useful information in understanding the 
hadronic effects on K* mesons. Therefore, we also take 
into account interactions of the kaon with light mesons 
during the hadronic stage of heavy ion collisions. 

To this end, we introduce effective Lagrangians for 
interactions between light mesons. The effective La- 
grangian methods have been used to calculate the scat¬ 
tering cross sections between J/'i/; and hadrons in or¬ 
der to estimate the amount of J/iji suppression in the 
hadronic matter [T 5 l - [l 8 l| . Recently, similar approaches 
have been applied to investigate the hadronic effects on 
the abundance of exotic mesons such as Dsj(2317) (Toj l 
and A(3872) mesons [2^[2l|. 

This paper is organized as follows. In Sec. II, we first 
consider interactions of both the K* meson and kaon with 
light mesons. Then we evaluate the absorption cross sec¬ 
tions of both mesons in the hadronic medium using ef¬ 
fective Lagrangians. In Sec. Ill we investigate the time 
evolution of the K* meson abundance by solving the ki¬ 
netic equation. In Sec. IV, we argue the important roles 
of the abundance ratio of K* mesons to kaons in heavy 
ion collisions. Sec. IV is devoted to conclusions. 

We have used throughout the work the isospin aver¬ 
aged mass for all hadrons, based on experimentally mea¬ 
sured masses [l^, e.g., rriK = 495.645 MeV. 

II. HADRONIC EFFECTS ON A* AND K 
MESONS 


K* meson produced at the chemical freeze-out can be 
absorbed or even produced through interactions between 
mostly light mesons during the expansion of the hadronic 
matter. We consider here K* meson interacting with 
the pions, p, K, and A*; A*7r —>• pK, K*p —>• ttK, 
K*K —>■ pTT, K*K* —> TTTT, and K*K* —>■ pp. The di¬ 
agrams representing each process are shown in Fig. [T] 
We introduce the following Lagrangians to describe the 
interaction between the K* meson and other two mesons; 


^pKK 

C-pK^K” 


= • (Ai9^7r-i9^A7r)-kH.c., 

= igpKK{Kfd^K - d^Krk) ■ 

= igpK^K^ [{dpR^'^TKl - K*''Td^kz) ■ fT 
+ {K*-T-d^p,-d^K*-T-p,)K*^ 

+ K*^{f-(rd^kz-f-d^irK:% (1) 


obtained from free pseudoscalar and vector meson La¬ 
grangians by introducing the minimal substitution. In 
Eq. (dl), A = (A0,A+) and K* = (A*°,A*+) de¬ 
note strangeness pseudoscalar and vector meson dou¬ 
blets, respectively, and tt and p denote the pion and p me¬ 
son isospin triplets, respectively, with Pauli matrices r. 
gTTK*K, gpKK, and gpK*K* are strong coupling constants, 
for which we use the empirical values, g-j^K-K = 3.25, 
gpKK = 3.05 [l^. We apply the SU(3) flavor symmetry 
to obtain gpK*K* = g-^K*K = 3.25. 

Using the above interaction Lagrangians we evaluate 
the amplitudes for all processes shown in Fig. [TJ The 
amplitudes of the A* meson absorption by pions, p. A, 
and A* mesons, without isospin factors and before sum¬ 
ming and averaging over external spins, are represented 
by 

M K*7z^pK 
M K* p^ttK 

iC * A * —V TT TT 

kiK'K'^pp 


= 

= 


We first investigate hadronic interactions of a A* me¬ 
son during the hadronic stage of heavy ion collisions. The 


where the amplitudes for the first A*7r —>■ pA and the 
second process A*p —>• ttA are 


J 


A ylio-) a */: 

My, = g-„K*KgpK*K*^iy 


1 


t — m^. + irriK^^K 


:( 




{pi -p^Tipi -psY 


‘‘K- 


{P2 +P4)p 


X ((2pi - pYipgau - {Pl +P?,)ugaP - {Pl - 2p3)aS'/3i/) , 
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FIG. 1: Born diagrams for the K* meson absorption by pions, p, K, and K* mesons. K*tt —>■ pK, (a) and (b); K*p —>■ irK, 
(c) and (d); K*K —>■ pn, (e) and (f); K*K* —>■ tttt, (g) and (h); K*K* —>■ pp, (i) and (j). 


and 


= -g^rK^KgpKKeieT - + 2F'2)m(F’ 3 + 2p4)i., 


s — m 


K 


= -g^K*KgpKKf^^i2-. -r(Pl “ ‘^P3)p{‘^PA-P2)^ 

L TTZ^ 


= —gTTK*KgpK*K*f^l(-2- 


tX 




{Pl +P2T{P1 +P2y 


s — rn\, + imjcTR* 

X ((2pi + P2)pgay - {Pl - P2)ugal3 “ {Pl + 2p2)aff/3i/) , 
respectively. Similarly, amplitudes for processes K* K —> pir and K*K* —> tttt are 


•'K* 


■) (P3 - Pa)p 


( 3 ) 


( 4 ) 


and 


1 y(e) O' *,5 

My, = g^K'KgpK'K-eies 


t — m^. + imK*T K- 


■( 




{pi -psTipi-psY 


i-K* 


■)(P2 +P4)a 


X {pPZ - Pl)agpy - ipi +P3)i^gap + (2pi -p3)/3ffai/), 


= gi^K-KgpKK^iey 


u — m 


-{2p4-pi)y2p2-p3)^, 


K 


M^p, = glK-K^i^T—ITiPi - 2P3)m(P 2 - 2p4)^ 




KaW — „2. p V 

JvIr* — g-KK*K^1^2 




respectively. Finally, the amplitudes for K*K* —>■ pp are 


-ipi - 2 pa)p{p2 - 2p3)i/, 


'•if 


( 5 ) 


( 6 ) 


A4CO _ 2 a */3 7 *5 

— gpK*K*^1^3 ^2^4 


1 


t — rriR, + imK*^K* 




(Pl -P3Y{P1 -P3Y 


LR* 


X (i2p3 - Pl)agpp - (pi +P3)/rPaP + (2pi “ P3)pflai.) ((P2 + P4)7fl5!. + (P2 “ 2p4)i.P75 + (P4 “ 2 ^ 2 ) 5571 /), 

4.,. , (Pl -P4)''(P1 -P4)'' 


((f) _ „2 ^a^*P^7^*5_5^_/ 

r* V 


_ 2 a *p 1 *c 


It — TO^. + imK^TR, 


‘‘K* 


X ((2p4 - Pl)a5/3p - (Pl + Pi)pgaP + (2pi - P4)p5ai.) ((P2 + P3)jgSi^ + (P2 “ ‘2P3)ugjS + (P3 “ 2p2)sg'ySj -(J) 
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In the above equations, pi denotes the momentum of par¬ 
ticle i. We keep the convention that particles 1 and 2 
stand for initial-state mesons, and particles 3 and 4 final- 
state mesons on the left and right sides of the diagrams, 
respectively. The Mandelstam variables s = (pi +P 2 )'^, 
t = {pi — and u = {pi — p 4 )^ have also been used. 
We apply here the K* meson propagator with its decay 
width, Tif., and use the isospin averaged value for the 
K* meson decay width, Tk* = 49.1 MeV [2^. 

In order to take the finite size of the hadron into con¬ 
sideration when evaluating amplitudes, we apply the fol¬ 
lowing form factor at each interaction vertex for the u, t- 
channel and the s-channel, respectively. 




A2-: 


A2 


cp 


Fs{q) 


+ wL 

A2 -I- 0^2 


( 8 ) 


with (p being the squared three-momentum transfer for 
t and u channels, and uP being the total energy of the 
incoming particles for s channel taken in the center of 
mass frame, niex is the mass of the exchanged particle 
in each diagram shown in Fig. ^ We set the cutoff 
parameter A to be A = 1.8 GeV [^. The final isospin- 
and spin-averaged cross section is given by, 


1 \Pf\ 

647r2s5i52 |p*| 




(9) 


where gi and g 2 are the degeneracy factors of the ini¬ 
tial 1 and 2 particles; gi = (2/i -|- l)(25'i -|- 1) and 
g 2 = (2/2-I-1)(25'2-I-1), respectively. |A4|2 represents the 
squared amplitude of all processes in Eq. © obtained 
by summing over the isospins and spins of both the ini¬ 
tial and final particles. \pi\ and |p/| in Eq. © stand for 
the three-momenta of the initial and final particles in the 
center-of-mass frame. 

Using the same method, we investigate hadronic effects 
on a AT meson during the hadronic stage in heavy ion 
collisions. We consider interactions of the K meson with 
pions, p, if, and K* mesons; ATtt —>■ pK*, Kp ^ irK*, 
KK —>■ TTTT, KK —>■ pp, and KK* irp. Among these, 
however, two processes Ktt pK*, Kp —>■ ttAT*, are 
same as the inverse processes of the K* meson interacting 
with p mesons and pions as shown in Fig. [TJ (c) and (d), 
(a) and (b), respectively, and the process KK* —>• np is 
same as that of the K* meson interacting with K mesons, 
diagrams (e) and (f) shown in Fig. [T] Therefore, all we 
need to consider more are the following amplitudes, 

■P^KK^pp = (10) 

for processes KK ^ tttt and KK ^ pp. We show the 
diagrams for theses processes in Fig. O 

The amplitudes for processes KK ^ tttt and KK ^ 
pp are 




9ttK-K 


t — m\, + imK»T K* 


(Pl +P3)/.(P2 +P4)i 



FIG. 2: Born diagrams for the K meson absorption by K 
mesons: KK — >■ tttt, (a) and (b); KK — >■ pp, (c) and (d). 


X ( - 5^" + 


{pi - psTipi - psY 


‘if* 






and 


X ( - + 


M)k = . "2 -A3)m(2A2 -Pa)^, 

Z 


{pi +Pi)p{p2+P2),. 
(pi -P4y{pi-P4)''\ 


u — rn]^, + irriK*^ K 


m 


(11) 


K* 


— 9pKK^4: ^3 


u — m 


-(2pi - P4)p.{2p2 - P3)„ 


K 


( 12 ) 


respectively. Then, using Eqs. © and © we evaluate 
the K meson absorption cross sections. 

Lastly we consider the possibility of the K* meson for¬ 
mation from pions and kaons. The scattering cross sec¬ 
tion for the K* meson production is given by the spin- 
averaged relativistic Breit-Wigner cross section; 




9 k - 47r __ 

9Kg^ Plm (mK* - i/sY + ’ 


with p^, gx, and gx* being the degeneracy of pions, 
K, and K* mesons, pi = (25^ -I- 1)(2U + 1), respec¬ 
tively, and Pcm the momentum in the center-of-mass 
frame. is the total decay width for a reaction 

Ktt — ?► K* —>■ Ktt as a function of ^/s. We take the 
following Y^-dependent decay width of the K* 

meson; 


Tx>^Kt,{Vs) = ^^^^PcmiVs), (14) 

We show in Fig. [3]the cross sections for the absorption 
of both the K* meson and the K meson by pions, p, K, 
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1/2 1/2-^ 

s -s„ (GeV) 


FIG. 3: Absorption cross sections for (a) the K* meson by 
pions, p, K, and K* mesons via processes K*-k —>■ pK, K*p —>■ 
tvK, K*K pvr, K*K* -/■ tttt, and A*A* -/■ pp, and those 
for (b) the K meson via processes KK — >■ tttt, KK — >■ pp, and 
Ktv K*. 


processes considered here, K*t: —>■ pK and KK ^ pp 
show very broad peaks above the threshold energy. 

We notice that K* meson is absorbed more by pions 
than by p mesons; the absorption cross section of the 
K* meson by pions, K*tt —>■ pK is larger than that by 
p mesons K*p —?► ttK. We also see that the annihilation 
cross sections for both the K* meson and K meson are 
larger when p mesons are produced than when pions are 
produced; the cross section for K*K* —>■ pp is larger than 
that for K*K* —>■ tttt, and the cross section for KK — pp 
is also larger than that for KK ^ tttt. 

The cross section for K*K* —pp is an order of mag¬ 
nitude larger than other processes, and seems to reflect 
the effect from two interaction mechanisms between three 
vector mesons. All particles participating in the process 
K*K* pp are vector mesons, and thus two C^k*k* in 
Eq. ([T]) are needed to describe the process K*K* —>• pp. 
It has already been shown that the interaction between 
three vector mesons increases the absorption cross sec¬ 
tion in the effective Lagrangian approach [i3- 

It seems awkward to observe that the cross section 
for K*K* —> pp should rise with increasing energy even 
though the form factor has been correctly used to kill 
the artihcial growth of the cross section with the energy. 
This behavior reminds the rise of the total cross section 
for pp collisions at high energy. It has been already well 
known that the resonance exchange is largely responsible 
for an increase of the cross section in high energy scatter¬ 
ing. In this study the K* meson exchange in the reaction 
K*K* —>■ pp causes the rise of cross section even at rel¬ 
atively low energy less than 1 GeV. The introduction of 
the decay width in the propagator Tif., however, does 
not contribute to this behavior at all. Instead it merely 
reduces a little bit the amplitudes for process having a 
K* meson exchange. Finally, we also find that the cross 
section for the formation of the K* meson from pions 
and K mesons, Eq. m is not small at all, compared to 
cross sections for other processes. 


and K* mesons via processes shown in Fig. [T]and Fig. [2] 
as functions of the total center-of-mass energy above 
the threshold energy Sq of each process. We see in Fig. 
[3] the general pattern that the cross sections have a peak 
near the threshold energy for the endothermic processes, 
e.g., K*tt —>• pK and KK —>• pp, while the cross sections 
for the exothermic processes, e.g., K*p —>• ttK, K*K —>• 
pn, K*K* — ^ TTTT, K*K* — >■ pp, and KK — ^ tttt, become 
infinite near the threshold. However, two endothermic 


III. TIME EVOLUTIONS OF THE K* AND K 
MESON ABUNDANCES 

We consider the time evolutions of the abundance for 
both the K* meson and kaon based on the cross sections 
evaluated in the previous section. We build a coupled 
evolution equation for both particles consisting of densi¬ 
ties and abundances for mesons participating in all pro¬ 
cesses shown in Fig. [TJ pions, p, K, and K* mesons. 


diVif.(r) 


dr 


= {crKp^K»7rVKp)np{T)NK{T) - {aK*-K^KpVK*-K)nT,{T)NK* {t) 

+ {aK7v^K*pVK^)n^{T)NKiT) - {crK*p^K7vVK*p)np{T)NK*iT) 
+ {<^p 7 v^K-KVp^)n^{T)Np{T) - W.if-^p^r//c.i?)niy(r)iViy.(r) 
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{^TTTZ^K* K*'^'7T7T )n,r('r)iV7r('r) — {ax*K*^TVTT'‘^K*K*)^K*{'’')^K»{T) 

+ {^pp^K-K-Vpp)np{T)Np{T) - {crK-K-^ppVK-K-)nK-iT)NK-{T) 

+ {a^K^K-v-^K)n-^ir)NK{T) - {TK-)NK*iT), 

= - {<JKK^^7vVKK)nK{T)NK(,T) 

+ {<^pp^KKVpp)np{T)Np{T) - {(JKK^ppVKK)nK{T)NK{T) 

+ (cr_R-7r-i._ftrpWif.^)n^(T)iV_R-(T) - {aKp^K--KVKp)np{T)NK{T) 

+ {<JK-p^K-^VK-p)np{T)NK>{T) - {(TKtt^K- pVK-K)nT,{T)NK{T) 

+ {<^p7v^K-KVp^)n^{T)Np{T) - {0K»K^p7rVK*K)nK(T)NK-{T) 

+ {Tk*)Nk*{t) - {a^K^K*v.„K)n^{T)NK{T), (15) 


where ni(r) is the density of a light meson i in the 
hadronic matter at proper time r, and Njir) is the abun¬ 
dance of the other light meson j in each process shown 
in Fig. [T] at proper time r. ni{T) for pions and p mesons 
is evaluated from 



by assuming that they are in thermal equilibrium, and 
varies in time through the temperature profile introduced 
below, Eq. m- We obtain Nj^r) by multiplying Eq. 
(fTBll by the hadronization volume V(t). In Eq. (ITSl) . 
Pi is the degeneracy factor for a particle i and K 2 the 
modified Bessel function of the second kind. 

ni(T) and fVj(r) are functions of the proper time 
through the temperature profile developed to describe 
the dynamics of relativistic heavy ion collisions. We 
use the schematic model of a system with an acceler¬ 
ated transverse expansion based on the boost invariant 
Bjorken picture 

V (t) = 7r[i?c + Vc{t - Tc) + ac/2(T - Tcf]^TC, 

T(T}=T,-(n-Tf)(^^^) , (17) 

\Tf -ThJ 

with T}i and Tf being the hadronization temperature and 
the freeze-out time, respectively. Eq. dnl) describes 
the system of the quark-gluon plasma expanding with 
its transverse velocity Vc and transverse acceleration Oc 
starting from its final transverse size Rc at the chemi¬ 
cal freeze-out time The temperature of the system 
decreases from the hadronization temperature to the ki¬ 
netic freeze-out temperature Tf. The values used in Eq. 
(1171) are summarized in Table HI 
In the rate equations, Eq. (fT51) . {oab^cdVab) is the 
thermally averaged cross section for initial two particles 
in a two-body process ab —?► cd given by 

i^ab—^cd'^ab) 


TABLE I: Values for the volume and temperature profiles in 
the schematic model Eq. Ha- 


Temp.(MeV) Time (fm/c) 
Rc = 8.0 fm Tc = 175 Tc = 5.0 

Vc = 0.4c Th = 175 Th = 7.5 

ac = 0.02cVfm Tf = 125 t/ = 17.3 


_ 1 / d^Pad^PbfaiPa)fb{pb)o'ab^cdVab 

l + 6ab J d^Pad^Pbfa{Pa)fb{Pb) 

1 

1 -I- 5ab ^rnlK2{ma/T)rrilK2{rnb/T) 

pOO 

X / dzKi{z)a{z'^T'^) 

J Zo 

X - {ma -f TOb)^/T^] - (toq - TOb)^/T^] ,(18) 

with Zo = max((ma +mb)/T, (rric + rad)/T), Ki and K 2 
being the modified Bessel function of the second kind, 
fi being the Boltzmann momentum distribution of the 
particle i, fi{p) = e“V^+^, respectively. Vab is the 
relative velocity of interacting particles a and b, Vab = 
\J{Pa ■ PbY - tnlmll{EaEb). (Tk*) in Eq. dHI) is the 
thermally averaged decay width of K* mesons, Eq. da, 
(r^.) = TK>{mK*)Ki{mK*/T)/K 2 {mK»/T), which has 
been obtained in the same methods as used in Eq. da- 
The K* meson abundance at t, Nk* , depends not only 
on the dissociation reactions like K*tt pK, K*p —>■ 
7rif, K*K —>■ pTT, K*K* —^ TTTT, and K*K* —>■ pp, but 
also on the production reactions, or the inverse reactions 
of the dissociation reactions, such as pK K*Tr, ttK —>■ 
K*p, pn K*K, TTTT ^ K*K*, and pp K*K*. We 
have taken both reactions into consideration in building 
the coupled equation for both the K* meson and kaon 
in Eq. da- We have used the detailed balance relation 
when evaluating thermally averaged cross sections of the 
inverse reactions from the results for forward processes 
shown in Fig. [3] The results are shown in Fig. U) 
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FIG. 4: Thermally averaged cross sections for the absorption of (a) a K* meson by pions, p, K, and K* mesons via processes 
K*-k —>■ pK, K*p —>■ -kK, K*K —>■ pTv, K*K* —>■ tttt, and K*K* —>■ pp, and those for (b) a K meson via processes KK —>■ tttt, 
KK —>■ pp, and Kn —>■ K*. Thermally averaged cross sections for their inverse processes pK —>■ K^tt, tvK —>■ K*p, pir —>■ K*K, 
TTTT —>■ K*K*, and pp —>■ K*K* for (c) a K* meson, tttt —>■ KK and pp —>■ KK for (d) a K meson. 


As we see in Fig. 21 thermally averaged cross sec¬ 
tions of the dissociation reactions are bigger than those 
of the production reactions for the exothermic reactions. 
In the case of the endothermic reaction like Ar*7r —>■ pK, 
thermalized production cross section is bigger than that 
for dissociation reaction. Both themalized cross sec¬ 
tions are comparable for the other endothermic reaction, 
KK pp. 

In general, the smaller threshold energy, mass, and de¬ 
generacy, the bigger thermally averaged cross section in 
the two-body process. We find that the thermally aver¬ 
aged cross section for K* formation, Ktt —>■ K* becomes 
more significant than those for other reactions. However, 
the unusually rising cross section in energy for the reac¬ 
tion K*K* —>■ pp has been suppressed in the thermalized 
medium as shown in Fig. Sla). 

When solving the coupled differential equation for both 
the K* meson and kaon abundances, we have treated 


abundances of their antiparticles K* and K mesons also 
as variables using the strangeness chemical potential p.^', 
i.e., Nk* and same for antikaons. In 

other words, we have not considered that K* mesons and 
kaons are in thermal equilibrium during the expansion of 
the hadronic matter, while we calculate the thermally 
averaged cross section Eq. (fT^ using the thermal distri¬ 
butions of hadrons involved. However, the initial yield 
of kaons at chemical freeze-out has been evaluated to be 
88.1 using the statistical hadronization model, Eq. (HD 
with the strangeness chemical potential ps = 10 MeV and 
the hadronization volume Vh = 1908 fm^ [2^, whereas 
the initial yield of K* mesons has been obtained from. 




dm 


mth ^BW {m - 
p'^dp 


Tr- 

mK-)^+Tl,/4: 

(19) 


0 _ 2 














FIG. 5: Time evolution of the K* meson and kaon abun¬ 
dances during the hadronic stage in central Au-|-Au collisions 
at y/SNN = 200 GeV. Ratio of the K* meson abundance to 
the sum of the K* meson and kaon abundances is shown in 
the inset. 


to take the width of the K* meson into consideration. In 
Eq. (fTOl) . rrith is the threshold energy for K* —>• Ktt de¬ 
cay channel and Nbw is the normalization constant for 
the Breit-Wigner distribution. We obtain the K* meson 
initial yield to be 55.7, which is slightly larger than 52.4 
obtained without including the K* meson width calcu¬ 
lated with the formula given in Eq. (HU). 

In Fig. [5l we show the abundances of the K* meson 
and kaon as a function of the proper time during the 
hadronic stage of heavy ion collisions at ^/snn = 200 
GeV. As we have expected, the K* meson abundance 
decreases due to both interactions of K* mesons with 
other hadrons and the decay of the K* meson to the pion 
and kaon, eventually becoming 35.6 at 9.8 fm/s after the 
chemical freeze-out. On the other hand the abundance 
of the kaon increases up to 106.5 at the end of hadronic 
expansion. We find that throughout the time evolution 
the sum of the K* meson and kaon abundances changes 
slightly from 143.8 to 142.1. We also show in the inset of 
Fig. [5] the variation of the ratio of the K* meson abun¬ 
dance to the sum of the K* meson and kaon abundances. 
The ratio decreases from the initial ratio from the statis¬ 
tical hadronization model 0.39 to 0.25 in the end. 

Based on the analysis we find that about 36% of K* 
mesons produced at chemical freeze-out disappear dur¬ 
ing the hadronic stage in heavy ion collisions, making 
the invariant mass reconstruction of the total K* meson 
difficult. We further find that the hadronic interactions 
shown in both Fig. [T] and Fig. [2] explain about 6% of 
the K* meson loss, and the K* meson decay is largely 
responsible for the K* meson reduction in the hadronic 
medium. Our result is comparable to the 30% reduction 
of the previous statistical model prediction 0.33 ± 0.01 



FIG. 6: A comparison between the K* meson and kaon abun¬ 
dances due to all hadronic interactions shown in Fig. [T]and 
Fig. [2] and those from the thermal model prediction at each 
time and temperature. 

to the experimental measurements 0.23 ± 0.05 p^ . 

We also consider the possibility of both the K* meson 
and kaon thermalization during the hadronic expansion. 
Assuming that both mesons are in thermal equilibrium 
with the hadronic medium we evaluate the K* meson 
and kaon abundances in time using Eqs. (HU), (HSl), and 
H3, and show the results represented by and 
in Fig. [6] 

As we see, keeps decreasing all the time. However, 
increases at the beginning of the hadronic stage, 
and finally decreases. This is due to the competition 
between the volume expansion and the decreasing rate 
caused by the thermal effects in Eq. (HSU, through the 
factor mif/T(r) inside the modihed Bessel function of 
the second kind, implying that and depends on 
the size and also on the lifetime of the expanding fire¬ 
ball, Eq. (IT71) . Nevertheless, one should note that the 
volume of the system expands in time with the total en¬ 
tropy almost preserved [^. We find that abundances of 
most hadrons decrease during the hadronic expansion in 
the statistical hadronization model but that of the pion, 
the lightest hadron, increases in time to compensate the 
loss of entropy from heavier hadrons. We argue that the 
same mechanism is happening to strangeness hadrons for 
some time during the hadronic expansion. However, the 
ratio -|- ) is not affected by the volume, 

and it keeps decreasing from 0.39 to 0.22 at the kinetic 
freeze-out. 

We notice that recent measurements of the K* meson 
yield in Pb-|-Pb collisions at Vsy y = 2.76 TeV at the 
Large Hadron Collider (LHC) provide 0.19 ± 0.05 as 
the /K~ ratio. This value is also smaller than the 

statistical hadronization model prediction 0.30 evaluated 
with the hadronization temperature 156 MeV [2^ at the 
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LHC energy. The measurements indicate that more K* 
mesons are lost during the hadronic expansion at LHC, 
leading to 37% reduction of the ratio. 


IV. THE ABUNDANCE RATIO OF K* MESONS 
TO KAONS IN HEAVY ION COLLISIONS 

Since the interactions of K* mesons and kaons with 
light mesons considered in Fig. [T]and Fig. [2] take place 
during the hadronic stage at both RHIC and LHC, it is 
necessary to understand general features of the variation 
of the K* meson abundance in heavy ion collisions. In 
order to analyze the reduction of the K* meson in the 
hadronic medium we simplify the coupled equation, Eq. 
(II3) by keeping the linear terms in Nk and Nk* only. 



dNK ( t ) ^ - jK-NK-ir), 

dr 

= -jkNk (t) + 7 if. Nk* (t) , ( 20 ) 

dr 

with 


IK* = {<^K*p^KTrVK*p)np + {(JK* 7 r^KpVK* 7 r)nTr 

+ (ric*), 

iK = {(TKTr^K*pVK-K)nTr + {(T K p^K*- kV K p)n p 

+ (o’iC7i-^a'*WiC7i-)?^7r- ( 21 ) 

When the thermal cross sections and densities of light 
mesons are independent of time, following analytic solu¬ 
tions are obtained from the coupled equation, Eq. (|20l). 

Nk*{t) = —N° + (n^, - 

7 V 7 / 

Nk{t) = 

'y \ y / 

where the initial yields for both hadrons, and N^, 
have been assumed, and iV° is the sum of the K* meson 
and kaon yields, iV° = -\- TV)*., at chemical freeze- 

out. The 7 in Eq. (l22ll is the sum of the K* meson 
and kaon widths in the hadronic phase, 7 = *)k + Ik* ; 
7 if • and jk play roles of the collisional broadening of the 
width of the K* meson and kaon in the hadronic medium, 
respectively. 

In Eq. time-independent terms represent abun¬ 

dances when time goes to infinity, and the sum of two 
solutions is preserved as its initial value . As 

time goes on Nk increases while Nk* decreases, and the 
rate at which the final number is reached in Eq. (l22l) is 
determined by the 7 which takes into account hadronic 
interactions of K* mesons and kaons with light mesons. 
If the 7 is large, the abundance can change significantly 
for a short time. 

Let us now investigate the time evolution of yield ratio 
of K* mesons to kaons from the analytic solution of Eq. 



FIG. 7: (a) Variations of *)k, *)k*, and 7 = 7 k -I- 7k* in time 
during the hadronic stage, (b) A comparison of the abundance 
ratio variation Nk*/{Nk + Nk*) obtained numerically from 
Eq. (fTSll for RHIC and R{t — 00) = 7k/7 from Eq. (EHl) 
evaluated at each time and temperature. We show in the 
inset how 7 (r — r^) changes in time. 


m, R{t) = NK*(.T)/{N*KiT) + NKir)); 


_ Nk* (r) _ Nk* (t) 

Nk* (t) + Nk (t) N° 

_ I / ^K* _ 'yK\ -*f(T-Th) 

7 V Afo 7 / 


(23) 


We notice that i?(r) is also composed of two parts; 
a time-independent part and a transient part. After a 
long time r the time-independent part i?(oo) = 7 k '/7 
is expected to represent the K* meson to the kaon ra¬ 
tio. How fast the abundance ratio approaches the time- 
independent part relies on 7 , the sum of the K* meson 
and kaon widths in the exponential function. 

With these in mind let us investigate the variation in 
the abundance of K* mesons and kaons obtained from 
Eq. (1151) . Since all thermally averaged cross sections and 
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densities of the light mesons in and ^k* are func¬ 
tions of a time, solutions of Eq. (ESI are different from 
the analytic solution of the simplified equation, Eq. (EOl). 
Nevertheless, we find that the solution of Eq. (HU keeps 
the same important characteristics of the analytic solu¬ 
tions from Eq. (|20l) . 

We first show in Fig. [71[a) 7 ’s obtained at Eq. 
m as a function of time. The 7 decreases in time 
from 0.33 c/fm to 0.26 c/fm as the system cools 
down from 175 MeV at th to 125 MeV at r/, re¬ 
flecting that interactions between hadrons become less 
vigorous as the temperature of the system decreases. 
The 7ic 5 K7t')^'k T (,^Kp^k*tt^Kp')^p T 

{(TKT,^K-VK7,)nT, in Eq. (EH decreases as the tem¬ 
perature decreases, but 'Jk* is almost constant dur¬ 
ing the hadronic stage; with decreasing temperature of 
the system, the part of 7 _ff, or {(Jk*p^k-kVk* p)np -I- 
{c^K*TT^KpVK*Tr)'n'Tr decreases whereas the thermal width 
of the K* meson (Tif.) slightly increases due to the fac¬ 
tor Ki{mK»/T)/K 2 {'mK*/T), meaning that K* mesons 
live shorter at lower temperature. 

We compare in Fig. [Hb) the abundance ratio variation 
/ {Nk + ) evaluated numerically from Eq. (fTSl) 

for RHIC to Jk/'J obtained from Eq. ( 1 ^ at each time 
and temperature, jk /7 represents the expected hadronic 
interaction width ratio between kaons and K* mesons 
plus kaons at each temperature and time. We anticipate 
that the abundance ratio of K* mesons and kaons in 
Fig. [7] approaches to as time goes on, like the 

ratio between those mesons obtained from the simplified 
rate equation, Eq. (1751) . We show in the inset of Fig. 
[7]Jb) how 7 (t — Th) varies in time at RHIC. As we see, 
7 (t — Th) increases up to 2.5 for 9.8 fm/c. We expect 
that the similar term with in the real solution 

suppresses the contribution of the time-dependent term 
as time goes on. The discrepancy between the abundance 
ratio and the in Fig. [TJ^b) is attributable to both 

the contribution from non-linear terms included in Eq. 
(1151) such as K*K pir, (e) and (f); K*K* tttt, (g) 
and (h); K*K* pp, (i) and (j) shown in Fig. [T] and 
KK TTTT, (a) and (b); KK pp, (c) and (d) shown 
in Fig. [7J and the time delay required to reach thermal 
equilibrium from the interactions of K* mesons with light 
mesons in the hadronic medium. 

Based on the above analysis we argue that the final 
ratio of the yield between K* mesons and kaons in heavy 
ion collisions is largely dependent on their interactions 
with other hadrons in the hadronic medium, jk and 
'jK* ■ Since 'y{T—Th) keeps increasing during the hadronic 
stage, the transient term with in the yield ra¬ 

tio R{t) plays a negligible role at a later time during the 
hadronic interaction stage. Therefore, we see that the 
relative interaction ratio jk /7 mainly determines the fi¬ 
nal yield ratio between K* mesons and kaons at the end 
of the hadronic stage. 

We show in Fig. [5] 7^/7 as a function of the temper¬ 
ature of the system. We also show in Fig. [5] measure¬ 
ments of the abundance ratio between K* mesons and 



FIG. 8: A variation of the 7 k/7 in temperature during 
the hadronic stage. We show measurements of the abun¬ 
dance ratio between K* mesons and kaons, Nk* /{Nk+Nk* ) 
0.23±0.05 at RHIC [3 and 0.19±0.05 at LHC gj 

kaons, Nk-/{Nk + Nk-), 0.23±0.05 at RHIC [3 and 
0.19±0.05 at LHC [T7| - We notice from Fig. [5] that the 
ratio of the K* meson and kaon in heavy ion collisions 
seems to reflect the interaction ratio between strange and 
light mesons, jk/j at the kinetic freeze-out temperature. 
We infer that the lower ratio A^k* /{Nk + Nk*) at LHC 
compared to that at RHIC is due to a lower kinetic freeze- 
out temperature at LHC than at RHIC. 

It has been argued that the degree of the reduction 
of K* meson yield during the hadronic stage in heavy 
ion collisions is attributable to a lifetime of the hadronic 
stage. Since the system of quark-gluon plasma at LHC 
is much larger than that at RHIC, it has been assumed 
that the lifespan of the hadronic stage at LHC is also 
longer compared to that at RHIC, and thereby more K* 
meson are lost in the hadronic medium at LHC. 

We find, however, from the investigation of the vari¬ 
ation in the yield ratio between K* mesons and kaons 
based on the solution of the rate equation, that the re¬ 
duction of the K* meson in heavy ion collisions reflects 
the interaction of K* mesons and kaons with light mesons 
at kinetic freeze-out. We argue that the degree of the K* 
meson abundance reduction in heavy ion collisions, or the 
reduction of the yield ratio between the K* meson and 
kaon, is largely attributable to the kinetic freeze-out tem¬ 
perature via the interaction of K* mesons and kaons with 
light mesons in the hadronic medium. The long lifespan 
of the hadronic stage just suppresses more a transient 
term, such as the second term in Eq. ESI, contributing 
little to the change of the K* meson to kaon ratio. 

As has been already mentioned, how fast meson abun¬ 
dances change in the hadronic medium is governed by 
the sum of all the interactions involved, i.e., the width 7 
in Eq. (1771) . Therefore, in addition to the hadronic inter- 
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actions considered in Figs. [T] and [2] all other interactions 
with various hadrons, i.e., nucleons, have to be taken 
into account to thoroughly understand the reduction of 
K* mesons in heavy ion collisions. Moreover, we also 
have to include more the feed down effects from heavier 
strangeness hadrons to fully consider the abundance ra¬ 
tio of the K* meson to the kaon, which are left for the 
future work. 


V. CONCLUSION 

We have studied the reduction of K* mesons in heavy 
ion collisions. We have focused on the hadronic effects on 
the K* meson and kaon abundances during the hadronic 
stage of the cental Au-|-Au collisions at ^snn = 200 GeV 
in order to understand the K* meson yield difference be¬ 
tween the experimental measurement and the statistical 
hadronization model prediction. We have evaluated ab¬ 
sorption cross sections for both kaons and K* mesons by 
pions, p, K, and K* mesons inside the hadronic medium. 
In describing the interaction between K* mesons and 
kaons and light mesons, we have introduced one meson 
exchange model with the effective Lagrangian. Further¬ 
more, we have built the coupled differential equation for 
K* mesons and kaons, and have solved it to investigate 
the time evolution of the K* meson and kaon abundances 
during the expansion of the hadronic matter. 

We have found that the K* meson and kaon abun¬ 
dances during the hadronic stage of heavy ion collisions 
are dependent on absorption cross sections and their ther¬ 
mal average. We have shown that the sum of K* and 
kaon abundances are almost preserved during the expan¬ 
sion, and the interaction of K* mesons with light mesons 
controls the reduction or production of K* mesons and 
kaons in the hadronic matter. Our analysis indicates that 
36% of the total K* mesons produced at the chemical 
freeze-out are lost during the hadronic expansion in the 


central Au-|-Au collisions at ^^snn = 200 GeV. We have 
found that among 36% about 6% of the total K* mesons 
are converted into light mesons by hadronic interactions, 
and the remaining 30% reduction is due to the decay of 
K* mesons to kaons and pions. We see that the loss of the 
K* meson abundance in the hadronic medium explains 
very well the discrepancy of the K/K* ratio between the 
statistical hadronization (lOl| m odel prediction and the 
experimental measurements [l3j. 

We have shown that the results obtained here can be 
applied to the analysis of the K* meson production at the 
LHC. We have found that all the interactions involved 
at RHIC must be present at LHC, and therefore widths 
7 if and jk evaluated at the RHIC energy can be ap¬ 
plied to the case at the LHC energy. Moreover, we have 
realized that the smaller ratio of K*/K measured at the 
LHC energy indicates a lower temperature of the kinetic 
freeze-out at LHC compared to that at RHIC. We have 
shown that the yield ratio between K* mesons and kaons 
is not mainly dependent on the lifetime of the hadronic 
stage in heavy ion collisions, and the hadronic interaction 
width ratio of strange mesons with light mesons, Jk/”/ 
determines the final yield ratio between K* mesons and 
kaons. We therefore conclude that studying the yield of 
the K* meson and its variation during the hadronic stage 
in relativistic heavy ion collisions provides a chance to 
understand not only the production of K* mesons but 
also the evolution of the hadronic medium in heavy ion 
collisions. 
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